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We consider the critical properties of the Bose-Fermi Kondo model, which describes a local moment 
simultaneously coupled to a conduction electron band and a fluctuating magnetic field, i.e., a dis- 
sipative bath of vector bosons. We carry out an e— expansion to higher than linear orders. (Here e 
is defined in terms of the power-law exponent of the bosonic-bath spectral function.) An unstable 
fixed point is identified not only in the spin-isotropic case but also in the presence of anisotropy. 
It marks the point where the weight of the Kondo resonance has just gone to zero, and the local 
moment fluctuations are critical. The exponent for the local spin susceptibility at this critical point 
is found to be equal to e in all cases. Our results imply that a quantum phase transition of the 
"locally critical" type is a robust microscopic solution to Kondo lattices. 

PACS numbers: 75.20.Hr, 71.10.Hf, 71.27. +a, 71.28.+d 



I. INTRODUCTION 

The interest in quantum criticality in strongly corre- 
lated metals arises primarily because it provides a mech- 
anism for non-Fermi liquid behavior. While the prob- 
lem appears to be important for a number of correlated 
electron systems - including high temperature super- 
conductors - the issues are particularly well-defined in 
heavy fermion metals Here, 
many materials have been shown to explicitly display 
a magnetic quantum critical point (QCP), including 
a growing list of stoichiometric or nearly stoichiomet- 
ric ones . A particularly striking puzzle has 
emerged from inelastic neutron-scattering experiments 
[ p~3| , p~4|JT5yT^ ] and magnetization measurements P, p^| , ^4| 
in some of these (nearly) stoichiometric materials. The 
dynamical spin susceptibility displays an uj/T scaling 
with a fractional energy/temperature exponent. In ad- 
dition, the same exponent is seen not only at the an- 
tiferromagnetic ordering wavevector but also essentially 
everywhere else in the Brillouin zone. 

In a recent work 1 17 T^|l^,^(| , a "locally-critical point" 
is found in Kondo lattice systems. This picture appears 
to explain the salient features of the aforementioned ex- 
periments. At a locally critical point, spatially extended 
critical fluctuations co-exist with spatially local ones. 
More microscopically, the criticality of the local Kondo 
physics is "embedded" into that associated with the mag- 
netic ordering of the lattice: the divergence of the spatial 
correlation length is accompanied by the destruction of 
the Kondo resonance. The microscopic analysis [ p7|Jl8| 
was carried out within an extended dynamical mean field 
theory (EDMFT), in which the Kondo lattice system is 
treated in terms of a single-impurity Bose-Fermi Kondo 
model supplemented by a self-consistency condition. The 
impurity model describes a local moment coupled at once 
to a conduction electron band and a dissipative bath of 
vector bosons; the bosonic bath describes a fluctuating 
magnetic field generated by the neighboring local mo- 



ments. The locally critical point is identified as a self- 
consistent solution when the Bose-Fermi Kondo model 
is treated to the first order in e within an e— expansion 
renormalization group (RG) procedure. Here e = 1 — 7, 
where 7 is the power-law exponent of the spectrum of the 
dissipative bosonic bath (defined in Eq. (|J)). 

There are several important questions that remain. 
First, does a self-consistent locally-critical solution arise 
at higher orders in e? Second, what happens to the lo- 
cally critical point in spin-anisotropic situations? Such 
anisotropy occurs in heavy fermions since the spin-orbit 
coupling is usually strong in these systems. 

These issues are addressed in the present paper. We 
carry out a detailed analysis of the Bose-Fermi Kondo 
model and show that an unstable fixed point exists not 
only in the spin-isotropic case but also in the presence of 
spin-anisotropy (both xy and Ising cases). In each case, 
this critical point describes a continuous transition from 
a Kondo phase, in which the local moment is quenched 
by the spins of the conduction electrons, to a "local- 
moment" phase where there is no Kondo resonance. For 
the isotropic and xy cases, we calculate the local spin sus- 
ceptibility to the order e 2 and, in addition, we determine 
the associated critical exponent, 77 (defined in Eq. (po|)), 
to all orders in e. In the Ising case, a slightly different RG 
approach turns out to be useful. |2l[] In all three cases, 
we find that the exponent r\ at the unstable fixed point 
is 



T} = e 



(1) 



This result turns out to guarantee the existence of a 
self-consistent solution of the locally-critical type in the 
Kondo lattice systems. 

We note in passing that the Bose-Fermi Kondo model 
is also of interest in other contexts. Historically, the 
model, in its Ising version, was first studied in the con- 
text of an EDMFT treatment of a spinless model in 
ref. 1 21 1, using an e— expansion. The spinful version was 
subsequently considered, also through an e— expansion, 
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in refs. [|22|,[23|. The Bose-only Kondo model was later 
extensively analyzed within a similar e— expansion in 
ref. p4 1 , which went to higher orders in e and established 
Eq. for a stable fixed point of that model to all or- 
ders. The Bose-only Kondo model appears implicitly in 
a single-impurity spin model first studied (using a large 
N approach) in ref. |2j| ; see also ref. |^6| . The Fermi-only 
Kondo model is of course a standard textbook problem 

The remainder of the paper is organized as follows. In 
Section we define the model and introduce the for- 
malism. Section [II is devoted to the RG analysis, to 



order e 2 , in the isotropic case. Section [V presents the 



explicit calculation of the local spin susceptibility to or- 
der e 2 , as well as the analysis of the associated critical 
exponent 77 to all orders. In Section we consider the 
effect of spin-anisotropy. Section [yj discusses the conse- 
quence of our results for t he lo cally critical solution in a 
Kondo lattice and Section VII provides a brief summary 
of our results. Some of the technical details are relegated 
to Appendices A, B C, and D. 



II. THE MODEL AND FORMALISM 

The Bose-Fermi Kondo model is defined as follows 
Wbfk = J S • s c + 2J E p cj, ff c pa 

p. a 

w P <l>p ■ 4>p- ( 2 ) 

p p 

A spin-i local moment, S, is coupled to both a fermionic 
bath (c p( j), through the Kondo interaction J, and a dis- 
sipative vector-bosonic bath (<p p ) with a coupling con- 
stant g. J is positive, i.e., antiferromagnetic, but g can 
be either positive or negative. (A sign change in g can 
be absorbed by a corresponding sign change in cf>.) The 
spectral function of the bosonic bath is taken to have a 
sublinear power-law dependence on energy, at sufficiently 
low energies: 

- w p ) - S(lu + w p )\ = (K%/n)\LupsgnLU (3) 

p 

for \u>\ < A. Here, the power-law exponent 

< 7 = 1 - e < 1. (4) 

The density of states of the conduction electron band 
near the Fermi energy is taken to be a constant: 



^2S(cj-E p )=N . 



(5) 



The fluctuating magnetic field competes against the 
Kondo-singlet formation. When the fluctuations are suf- 
ficiently slow such that e > 0, g is a relevant coupling in 



the RG sense and it can compete with the marginally- 
relevant J coupling. The physics of this competition is 
amenable to an e— expansion 22 2^j2l[]. This is in con- 
trast to the Kondo fixed point, which occurs at an infinite 
coupling and can only be addressed by strong-coupling 
methods f27|. 

Following Smith and Si I^Jl^Q, we adopt the 
Abrikosov representation of the spin in terms of pseudo- 
/-electrons [p8|, 



(6) 



where t x,v,z are the Pauli matrices. In this representa- 
tion, the Bose-Fermi Kondo Hamiltonian takes the fol- 
lowing form, 

n= Tin + Ttj + H. 



j-r n,g, 

Ho= fl U + ^Ep c\ a c pa + w p $J ■ (f p , 

p 



(7) 



ri.r- 



4 ^ 



GG 



where A is an energy level for the /—electron that will 
be set to 00 at the end, g = ±1, <j> = J2p($p + ft-p)} an d 

To analyze the critical behavior of the Bose-Fermi 
Kondo model, we carry out an RG procedure based on 
a dimensional regularization with a minimal subtraction 

(MS) of poles gl^H!!. We define a renormalized 
field / and a dimensionless coupling constant g by 



Sb = z) /2 f, 

9b = gZ^Zgfi^ 2 , 



(8) 
(9) 



where Zf is the wave-function renormalization factor for 
/ electrons, gs is the bare coupling constant, Z g is a 
coupling constant renormalization of g, and [L is a renor- 
malization energy scale. 

The "dimensional" regularization for the conduction 
electron density of states is done by introducing an e' 
through 



J2d(Lo-E p ) = N 



(10) 



Accordingly, the dimensionless coupling constant J is de- 
fined by 



J B = JZfZj/ 



(11) 



Note that e' is introduced strictly for the purpose of fa- 
cilitating the MS analysis; it is set to zero at the final 
stage of the calculation. 
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III. RG ANALYSIS 



In this section, we carry out the RG analysis to order 



A. RG flow equations 

It is known that, to the linear order of e, there is a crit- 
ical point at (K Q g*) 2 = e/2, (N J*) = e/2 H||,0@. 
To obtain the corrections to the order e 2 , we need to cal- 
culate the self-energy and vertex corrections beyond the 
orders J and g 2 and include also terms to the orders <? 4 , 
J 2 , and g 2 J. The details of our calculation are given in 
Appendix A. 

The renormalization factors are obtained as 



Z f = 



z„ 



Zr = 



15 

32c 2 
1 

+ 4e 
9 

321 2 

1 - -jNoJ 

e 

1 

+ 4e 



+ 



(Kog) 4 
(Kog) 2 
(Kog) 4 



8e 
1 

16e 
1 



-(Kog) 4 



(K g) 



(K g) 4 



^(NoJ? + ±- 2 (NoJ? 



32e ; 



■(K 09 y 



:(K g) 4 



(12) 



They can then be used to determine the beta functions 
for the coupling constants g and J, defined as |33| 



P(g) = ^\ aB j B , 

dJ 

P(J) = VJ^\ob,Jb- 

The results are given as follows 



P(g) = -9 ( ^ - (Kog) 2 + (Kog) 4 



(13) 



(NoJf 



(3(J) = -J (N J) 



(N Jf 



-j(-(K g) 2 + {K gf). 



B. The phase diagram 



(14) 



The RG equations (|l4|) yield several fixed points. 
There are two stable ones. The fixed point located at 



0.15 




0.05 



0.15 



FIG. 1. RG flows in g - J plane when e = 0.1. K (Kondo) 
and L (local-moment) denote the two stable fixed points. The 
separatrix is determined numerically. C denotes the unstable 
fixed point (the critical point). 

g* = and large J* ( "K" in Fig. |l|) , as usual, specifies the 
Kondo phase |p4| - Here, the local moment is quenched 
by the spins of the conduction electrons, leading to the 
development of a Kondo resonance. 

Another stable fixed point ( "L" in Fig. is located at 

1 



(N r) = 



(15) 



We will call it a "local- moment" fixed point, emphasizing 
the fact that it describes a phase in which no Kondo res- 
onance arises. The dynamics in the local-moment phase 
is controlled by the coupling of the local moment to the 
dissipative bosonic bath alone. 

A separatrix specifies the boundary of the domains 
of attractions for these two phases in the J-g parame- 
ter space. Lying on this separatrix is an unstable fixed 
point, or a critical point ("C" in Fig. |l|). It is located at 



(K g*) 2 = ^ + le 2 + O(e 3 ) 



(N J* 



) = 7; + 0(e 3 ) 



(16) 



The critical point marks the point where the spectral 
weight of the Kondo resonance goes to zero. It captures 
the competition between the Kondo coupling of the local 
moment to the conduction electrons on the one hand, 
and its coupling to the fluctuating magnetic field on the 
other. 

The RG flow can be solved numerically. In Fig. |[ we 
give the result for e = 0.1. 
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IV. THE LOCAL DYNAMICAL SPIN 
SUSCEPTIBILITY 

We now turn to the calculation of the local suscepti- 
bility. 



A. Critical susceptibility to order e 

Consider the local spin susceptibility to order e 2 . We 
first calculate the bare spin susceptibility, i.e., the au- 
tocorrelation function of the unrenormalized local spin 
operator specified by Eq. (||), 

Xb (t)^ 1(T t S-(t)S+(0)) = lim ~ XB (r), 



(17) 



For completeness, we also briefly discuss the local spin 
susceptibility at the local- moment fixed point. We as- 
sume that the susceptibility here also has a power-law 
form: 

xM-^(^)", (22) 

with an anomalous dimension r/L. 77^ has the same de- 
pendence on g*, J* as in Eq. (21). Even though the cou- 
pling constants g* , J* now take the values specified by 
Eq. (|l5|), it turns out that the e 2 terms again cancel with 
each other leaving 

r]L = e. (23) 
The critical amplitude is also given in Appendix B. 

B. Critical exponent to all orders in e 



where 5 ± = (S x ± iS v ). (For notational simplicity, 
we've dropped the subscript B in the operators used in 
Eq. @.) 

Correspondingly, we define a spin renormalization fac- 
tor Z which relates the renormalizcd spin susceptibility, 
X, to the bare susceptibility, 



(18) 



The calculational details are given in Appendix B. From 
Eq. (B6), Z is given by 



2(K gf , (K Q gf (A Q J) 2 
Z=l — + — ( 19 ) 

At the critical point, the local spin susceptibility is 
expected to have a power-law form: 



(20) 



for |t| ^ To, where To = 1/A is a cut-off scale. Here, 77 is 
the anomalous dimension. It can be calculated from the 
renormalization factor 



dlnZ, 



ls=9* ,J=J* 



dft 

n " ] — • 1U] —) >• 

2(Kog*) 2 ~2(K Q g*) 4 + (N J*) 2 
c. 



(21) 



It is striking that, although the order e 2 corrections are 
present both in the beta functions and in the critical cou- 
pling constants, they cancel out in the critical exponent 
77, leaving only a non-zero linear term. 

The critical amplitude is given in Appendix B. 



In this section, we show that 77 = e given in Eq. ( pl| ) 
is in fact exact to all orders. 

First, we explore the reason of the cancellation of the 
e 2 terms in 77 seen in Eq. ( pi] ) . Using the equations ( p^ ) , 
(p9|), we can easily verify that 

Z- 1 = {Zj l Z g f (24) 

both to linear and second orders in e. Combining Eq. 
(§J) with Eq. (g) lead to 

g B =gZ-^ l f'\ (25) 

Differentiating the logarithm of both sides with respect 
to In /j, keeping g B and J B fixed, we have 

, 2/%). 



n 



9 



\g=g*,J=J* 



At the fixed point, (3(g) = 0, which gives 
77 = e. 



(26) 



(27) 



There is also a consistency check for our result. We 
can start from Eq. (|ll]). Combining this equation with 
Eq. (H) lead to 



J B = JZ- 1 ' 2 (Z J /Z g ). 



(28) 



(We have set e' = 0.) Again, differentiating the logarithm 
of both sides with respect to ln/_t, and again keeping g B 
and J B fixed, we end up with 

dln(Zj/Z g ). 



V 



\g=g-,J=J* 



(29) 



To be consistent, the RHS should be equal to e. This is 
indeed verified to both linear and quadratic orders in e. 

We now proceed to higher order contributions. We 
show in Appendix B that Eq. (|24|) is valid to all orders. 
Eq. ( |26] ) is then valid to all orders in e. By extension, 
77 equals to e to all orders in e. Our reasoning basically 
parallels that of ref. \24\ for the stable fixed point in the 
Bose-only Kondo problem. 
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V. ANISOTROPIC CASES (XY AND ISING) 

So far, we have discussed the spin-isotropic case. We 
now turn to the effect of anisotropy in spin-space. 



A. RG equations 

We introduce separate parameters for the transverse 
and longitudinal spin couplings as follows: 

hj= J -f J2 ™'fUA>c„> + ^(fUi c h + H - c -)> 



n g = | £ vfiw + ^J\hr + h.c), 



(30) 



The RG procedure parallels that for the isotropic case, 
except that we need to differentiate between the longi- 
tudinal and transverse couplings in all the contributions, 
and associate different renormalization factors for the dif- 
ferent types of vertices: 



gf 
J? 



t/2 



(31) 



where i =_L, z. Up to the order of our interest, the results 
for the renormalization factors are given in Appendix C. 
The resulting RG equations are given as follows, 



(K g ± ) 2 + (K g z f 



(K g ± f(K 09z f + (K o 9± ) 4 
2 

{No J±) 2 + (N J Z ) 2 



P(j9z) = ~-j9z + {Kogi_) 2 g z 

- (K g^ 2 (K Q g z ) 2 g z 
^(NoJ±f 

H n 9z 



-9x 



(32) 



(33) 



(K og± ) 2 + (K g z f 



J± 



(K 09l _f{K g z f + (Kpg^ 
2 J± 

(N J z) J ± + ^lp^lj ± 



(3(J Z ) = (K Q g ± ) 2 J z 

- (K g±) 2 (K g z ) 2 J z 

- (No Jx) J± + &^J, 



(34) 




FIG. 2. Schematic RG flow diagram for anisotropic g 
coupling when J = 0. The Ising fixed-point and the xy 
fixed-point are stable while the isotropic fixed point is un- 
stable. 



Consider first the J = case, i.e. the Bosc-only Kondo 
model. The RG flow is given in Fig. g. There are two 
more fixed points in addition to the isotropic one. The 
xy fixed point is located at 



(K Q gl) 2 =e + e 2 +0(e 3 ) 
(Kog* z )= 0. 



(36) 



While the Ising fixed point is nominally located at g\ = 
and g* — oo. Both the isotropic fixed point and the xy 
fixed point are accessible by the e— expansion. The Ising 
fixed point, on the other hand, is beyond the reach of the 
perturbative RG scheme. 

The J— coupling can lead away from the "local- 
moment" fixed point not only in the isotropic case, but 
also in the xy and Ising cases. The isotropic case has 
already been discussed in the previous sections. We now 
turn to the xy and Ising cases, respectively. 



B. Critical behavior in the xy case 

We will now set 9z = in Eqs. @- ©. [Eq. @ 
implies that g z will stay at zero under the RG transfor- 
mation when its initial value vanishes.] An unstable fixed 
point occurs at 

(i^I) 2 =e+^ + 0(e 3 ) 

NoJ1= T2 + Tk/ + 0(eS) 



(35) 



n j:= 



0(e 3 



(37) 
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FIG. 3. Schematic RG flow diagram for the xy case. Here 
we show the projection of the RG flow in the three dimen- 
sional parameter space (J±, J z , g±) onto a) the Jj_ — g± 
plane and b) the J z — g± plane. 

Note that g z remains irrelevant near this fixed point, es- 
tablishing the consistency of our analysis. The RG flows 
are shown in Fig. [| 

The perturbative correction to the local spin suscep- 
tibility can be calculated as in the isotropic case. The 
most general expression for the susceptibility as well as 
its renormalization factor Z are given in Appendix C. 
The corresponding expression for 77 is also given there, in 



Eq. (C9). Using the values of the coupling constants at 
this fixed point, Eq. (|37|), we find that 77 remains to be e 
to order e 2 . 

The same argument made in Appendix B for the 
isotropic case carries through here for the xy-case as well, 
resulting in 



Z 1 — (Zf 1 Z g± ) 2 , 



(38) 



to infinite orders in e. So 77 = e is again valid to all orders 
in e. 

We now again briefly examine the local spin suscepti- 
bility at the local-moment fixed point in this case, where 
Jl = J* = 0, and 



(39) 



The anomalous dimension t]l has the same expression as 
in Eq. ( C12| ). Substituting the values of the stable fixed 
point, we find that r]L is equal to e, just as in the isotropic 



C. Critical behavior in the Ising case 

None of the non-trivial fixed points in the Ising case is 
within the reach of the perturbative RG scheme. For the 
unstable fixed point, for instance, setting g± = in Eqs. 

e that is still small, but 



will yield a (A^l) 2 
J* and g* that are of order unity. The latter violates our 
starting assumption. 




FIG. 4. Schematic RG flow diagram for the Ising case, pro- 
jected onto the plane spanned by yj — n g plane. 

This problem, however, has already been studied in 
rcf. A finite J z and g z can be handled by introduc- 
ing a so-called "kink-gas" representation. In this repre- 
sentation, the unstable fixed point is still accessible by an 
e— expansion. The calculation is outlined in Appendix D. 
The RG equations are no longer constructed in terms of 
the bare couplings. Instead, they are given in terms of 
the stiffness constants Kj , K g , and the fugacity yj , whose 
initial values are specified by the bare parameters as fol- 
lows, 



1 



K f=[l-^tan- 1 (^iV jf)] 2 



Hi 



r (7) 1- 



B\2 



(40) 



The RG flow projected onto the yj — n g plane is shown 
in Fig. |]. In the vicinity of the unstable fixed point, k.j 
is irrelevant. 

The exponent 77 at the unstable fixed point is calculated 
in Appendix D to the linear order in e. Again, it is found 
that 



i] = e 



(41) 



In the same Appendix, it is also shown that 77 is in fact 
the same as that for the critical point of a classical ferro- 
magnetic Ising chain with a long-range interaction that 
decays in distance in terms of a power-law exponent 2 — e. 
For this problem, it has long been held that Eq. ( f4l| ) is 
in fact exact ]3^ , |36[ |. 

In passing, we note that this analogy also allows us 
to state what happens to the local susceptibility at the 
local-moment fixed point shown in Fig. (Q). Here, x(r) 
picks up a constant (r— independent) piece. At the same 
time, the "connected" susceptibility (defined such that 



G 



the constant piece does not appear) decays with an ex- 
ponent rjL — 2 — e 



VI. LOCALLY CRITICAL POINT OF THE 
KONDO LATTICE 

We now discuss the consequences of our results for 
the locally critical point of the Kondo lattices. A self- 
consistent microscopic treatment of the Kondo lattice 
model has been presented in detail elsewhere |l7| , |l8| . 
This was carried out within the EDMFT approach de- 
veloped in pi| , p8| , p9[ , (This approach is a generalization 
of the dynamical mean field theory jl0| such that mag- 
netic fluctuations are also taken into account.) The self- 
consistent equation reads as follows, 



< x(q,w) > q = xiocM, 



(42) 



Here, the lattice dynamical spin susceptibility has the 
form 



x(q,w) 



1 



M{oj) + 7 q ! 



(43) 



where M(cu) is the "spin self-energy" and I q describes 
the exchange (RKKY) interaction between the local mo- 
ments. Eq. ([l2]) is simply a statement of translational 
invariance, as it equates the on-site spin susceptibility 
of the lattice system (LHS) to the susceptibility of any 
local moment (RHS). The locally critical solution arises 
when xioc(w) is singular - signaling the destruction of 
the Kondo resonance as discussed in previous sections 
- at the point where the peak susceptibility \(Qaf,^) 
just diverges. It is shown in rcfs. 0,^8| that one re- 
quirement for a locally critical point is that the magnetic 
fluctuations are two-dimensional. This can be easily seen 
from Eqs. ( f42||l3"| ) . The fact that the peak susceptibility 
x{Qaf,u) diverges means that M(0) = —Iq af - Recog- 
nizing that for q ~ Qaf, the interaction has the generic 
form, Iq = Iq af + a(q - Qaf) 2 , we have 



x(q~ Qaf 0) 



«(q - Qaf) 2 



(44) 



Now, the fact that the local susceptibility, xioc(<^ — > 0), 
is singular demands that the q— averaging of the LHS of 
Eq. (^) should diverge. In other words, integrating the 
RHS of Eq. (HJ) over q should yield a singular result. 
Two-dimensionality provides the phase space for this. 

This condition, however, is not sufficient. It turns out 
that a second condition must be met for the existence of 
such a locally-critical solution. The condition is precisely 
Eq. ([!]) . The origin of this second condition is somewhat 
subtle, but has to do with the fact that M(u>) and Xloej i n 
addition to being related through Eqs. (^,^), must also 
satisfy a Dyson-like equation of the (self-consistent) Bose- 
Fermi Kondo model. The details can be found in ref. fEs| . 



By proving that Eq. (nj) is valid to all orders in e, and 
also in spin-anisotropic cases, we have then established 
that the locally critical point is indeed a robust result 
within these microscopic considerations. (The arguments 
for the robustness of the locally critical point beyond the 
microscopic approaches can be found in refs. ]l7||lS| ].) 

We stress that a key element of the locally critical point 
of the Kondo lattice is that, the criticality of the local 
degrees of freedom is embedded in the criticality associ- 
ated with the long-range ordering in the lattice. In other 
words, the point where a local energy scale turns to zero 
coincides with where the spatial correlation length just 
diverges. This is very different from the self-consistent 
spin liquid solutions discussed in other contexts f25| , ^6f , 
which correspond to a phase instead of a critical point. 
This difference is also reflected in the fact that stabiliz- 
ing the locally critical solution in Kondo lattices requires 
two-dimensional magnetic fluctuations. 



VII. SUMMARY 

In short, we have carried out a detailed analysis of the 
Bose-Fermi Kondo model both when the spin-rotational 
invariance is satisfied (isotropic) and when it is broken 
(xy and Ising). In each case, we have identified a crit- 
ical point that separates a Kondo phase, where the lo- 
cal moment is quenched by the spins of the conduction 
electrons, and a local moment phase where there exists 
no Kondo resonance. This unstable fixed point marks 
the point where the spectral weight associated with the 
Kondo resonance has just gone down to zero. 

In all three cases, we find that the exponent for the lo- 
cal spin susceptibility at the unstable fixed point is equal 
to e. We note that the three cases correspond to three dif- 
ferent one-dimensional statistical mechanical problems. 
It is remarkable that, the susceptibility exponent at the 
unstable fixed point is insensitive to these differences. 
However, the same cannot be said about the stable fixed 
points: here, the susceptibility exponent for the Ising 
case is very different from the isotropic and xy cases. 

Our results have important consequences for the lo- 
cally critical behavior in Kondo lattices. In particular, it 
guarantees that a quantum phase transition of a locally- 
critical type [17 18| ll| is a robust microscopic solution 
to the Kondo lattices: it arises to all orders in e and, in 
addition, not only when the system is spin-rotationally 
invariant but also in spin-anisotropic situations. 
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APPENDIX A: RG PROCEDURE TO ORDER e J 

In this appendix, we give the details of the RG analysis. 
We adopt the dimensional regularization and minimal 
subtraction scheme [f|||],|3|,|3| . 

We first carry out a renormalized perturbation calcu- 
lation. The quantities of interest are the J— and g— 
vertices: 



Tj{uj + A)= Jfi e 7j(w + A) 



(Al) 



as well as the /—electron self-energy, + A). 

G-\lo + A) = u) - E/(w + A) = Z f Gg l {u) + A) (A2) 

Within the renormalized perturbation calculation, the 
bare parameters in Eq. (Q) are replaced with the corre- 
sponding renormalized ones. We end up with a Hamilto- 
nian containing one part that has the form of Eq. (g) (the 
tree level Hamiltonian) and another part which describes 
the counterterms: 

n ct = J2( z f - !) A ft fa + (Zj - i)j/s ■ S C 



-(Z g - l)g^ /2 S ■ 



(A3) 



There are then two types of perturbative diagrams, one 
coming from the tree level Hamiltonian, which we will call 
direct perturbative contributions below, and the other 
from the counterterms. 

The singular contributions are kept track of through 
poles as a function of e and e'. By demanding that such 
poles are "minimally" removed so that the renormalized 
quantities T g , Tj, and Gf are regular (as a function of e 
and e') , we fix the wavefunction and vertex renormaliza- 
tion factors Zf,Zj, and Z g order by order in the coupling 
constants. From these renormalization factors, we can 
determine the RG beta functions as well as the anoma- 
lous dimensions. Note that what we are carrying through 
is a double expansion in e and e'. 

In the following, we separate our discussions according 
to the order in J and g. 




(z f -ir 




FIG. 5. Self-energy diagrams for the /—electrons. The 
dashed, solid and wavy lines denote the propagators of the 
/-electrons, conduction electrons and vector bosons, respec- 
tively. The cross in (e) and dot in (f) denote the counterterm 
Zf — 1 and Z g — 1, respectively. 



1. First order result 

To the hrst order of perturbation, there are only cor- 
rections to the J— coupling from processes that involve 
the J— coupling alone. The diagrams are specified by 
Figs. 5(a) (b) in Ref. @]. The result is given by 



tS'V + a) 



No J. ii 



€' 

N J 



(— r 



N J\n( — ) 

— LU 



(A4) 



Adding the counterterm (Zj — 1) and demanding that 
the pole terms cancel, we find that 



{Zj 



1 



1 



-,{NoJ). 



(A5) 



2. Second order result 



To the order g 2 and J 2 , the contributions to the 
/—electron self-energy are shown in Figs. ||a and |^b and 
they yield the following expression: 



£f^ + A) = -c,f(^os) 2 ( — r 
' 4e — u) 



-u-?-(N jf(j^r', 

foe' — oj 



(A6) 



By demanding that these poles of the self-energy be can- 
celled by the counterterm contribution —(Zf — l)u>, we 
fix 
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(Z f -lp=-^(K 09 f~^(N J) 2 . (A7) 

To the same order, the corrections to the coupling con- 
stant g are specified by Fig. 6 of Ref. |j~8j and they yield 



1 



7^(c + A) = --(l- e )(X 5) 2 (^r 
" 4e —lu 

16e' — u; 

This fixes the second-order contributionn to Z„ 



1 



,(2) 



(A8) 



(A9) 



The corrections to the coupling constant J include not 
only the direct perturbation diagrams, but also contribu- 
tions from the first order counterterm (Zj — lY^J. Figs. 
5(c) (d) in Ref. G§ give 



1 

If 
1 

16e 



-U! 



7 (7V J) 2 (^) 2e '. 



(A10) 



The so-called "parquet" diagrams as in Kondo problems, 
shown in Figs. g(a)-(f), give 



t u —U) 



(AH) 



The diagrams with the counterterm (Zj — 1)^ J, shown 
in Figs. |(g)(h), give 



7 S 2)m (o; + A) = 2(Z J -l)« 7 «(a; + A) 



■( — ) e , 



(A12) 



where the prefactor 2 is a symmetry factor. From these 
contributions, we obtain the second order correction to 
Zj: 

(Zj 1) (2) = j- e (Ko9) 2 + Y^oJ) 2 + ^{N jf . 



3. Third order result 



(A13) 



To the third order of perturbation, we need only to con- 
sider the g 2 J order corrections to the J vertices. (There 
are also J 3 order corrections to the J vertices, but these 
contributions are higher than e 2 order.) The direct per- 
turbative diagrams are shown in Fig. ^. Summing up the 
contributions of these diagrams, we get 



iff- + "45 



e(e + e') 



x(X o5 ) 2 (7V J)(— ) e+e '- 

— UJ 



(A14) 





^ ,-7 « W 




FIG. 6. The vertex correction diagrams for the coupling J 
to the J 2 order. Here, (a)-(f) are the "parquet" diagrams; 
(g)(h) are diagrams with counterterms (Zj — l)' 1 '. 



Fig. H shows the diagrams with counterterms to g 2 J or- 
der. The contributions from Figs. ||(a)(b) give rise to 

1 f II {u J + \) = -(Z } -l)V 1 t l\u + \) 
3(K g) 2 (N J) . [i 



4ee' 



(A15) 



those from Figs. ||(c)(d) yield 

7 ( 3)/// (c + A) = 2(Zj - 1)^7^^ + A) 
_ (K g) 2 (N J). n 



2ee' 



(—r 

UJ 



(A16) 



and finally, Fig. ||(e) leads to 
7 ? )IV (« + A) = (Z J -1)W 7 ?)( W + A) 



h - h) (Ko 9 nN J)(^-y, (A17) 



Summing up j^p I ' II ' III ' IV an d subtracting poles, we 



have 



(^- 1 ) ,3> = -( ? (7^o + i)(^( JV « J >- 



(A18) 



Note that the simple poles in Eqs. (A14, A17) cancel with 
each other, leaving only double poles in the net result, 



Eq. (A18). Therefore, no g 2 J terms will appear in the 
beta functions. 



9 




FIG. 7. The direct perturbative vertex diagrams for the 
coupling J to the g 2 J order. 




FIG. 8. The vertex correction diagrams for the coupling J 
to the g 2 J order with counterterms (Zf — l)' 2 ', (Zj — lp 1 ' 
and (Zj - 1) (2) . 



7! if? ? 7? 



71 



'<z,-i>' ! 



5 ^ 



FIG. 9. The vertex correction diagrams for the coupling g 
to the g 4 order, (a)-(f) are the direct perturbative diagrams 
while (g), (h) are diagrams with counterterms 

4. Fourth order result 

We now turn to contributions of order g 4 , the only 
terms we are interested in to this order. First, we con- 
sider the /—electron self-energy. Figs. ||(c), (d) are 
the direct perturbative contributions to /—electron self- 
energy of this order. Figs. ||(e), (f), on the other hand, 
are the diagrams with counter terms of order g 2 and J 2 . 
Keeping only the g 4 order terms, we have: 



,(4)(| c +|d) 



(id + A) = LU 



4 4 W(u, + A) 



15 
32? 

9 



|)(iv 03 ) 4 (^ 



■(Z r l)( J )E?»(w + A) 



-to 



16e 2 

9 



(Ko 9 n—r 

— LU 



2{Z g - l)^Y,f{u + A) 



-UU^-{K,g)\^) 



(A19) 



■se' —lo 
The resulting corrections to Z ' f in this order is given by 



(Z f -l)^=-^(K g) 



15 



(A20) 



Next, we consider the g— vertex corrections. They are 
specified in Fig. 0. The contributions from the direct 
perturbative diagrams, Figs, ^(a)-(f), sum up to 



(A21) 



while those from the counterterms, Figs. |9|(g)(h), yield, 
7 W ^ico + X) = 3(Z g -l)^^\co + X) 
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FIG. 10. The vertex correction diagrams for the coupling 
J to the g 4 order. 




FIG. 11. The vertex correction diagrams for the coupling 
J to the g 4 order with counterterms. 



Therefore, the correction to Z g is given by 



(A23) 



Similarly, the J— coupling corrections contain both di- 
rect perturbative diagrams (Fig. [lO]) and diagrams with 
counterterms (Fig. [ll]) which are given by 



7, ( / )/ ^ + A)=(^- I ^) /« : — (A24, 
and 



^4/ M -j2e 
-U> 



lf"[u + A) = (2(Z g 1)( 2 ) - 2{Z, ir) 7 «( W + A) 

+ (Z, 7 - l)( 2 > 7 f ( W + A) 
9 



16e 2 



loe/ — 



respectively. The resulting correction to Zj is 
(Zj-l)W = ^(K g)*-j-(K gf. 



(A25) 



(A26) 



We can now collect the renormalization factors Zt, Z g , 
and Zj to the orders of our interest. The results are given 
in Eq. A 



5. The beta functions 

After obtaining the renormalization factors, we can cal- 
culate the beta functions, defined in Eq. (Ol). 



Taking the [i derivative of Eqs. (^"Tj) at the fixed bare 
couplings, we have 

8(a) dMZji-Zg) 



9 
8{J) 



On 



d\a{Z~ l Z g 



c 

2' 



dJ 

d\n{ZfZj) 



J ' dg 
dhx{Zj l Zj) 



d,J 



(A27) 



After some algebra, and setting e' = 0, we obtain the 
beta functions of g and J as in Eq. (jlj) . 

6. Renormalizability 

We close this appendix by commenting on two issues 
regarding renormalizability. Our perturbative results for 
the self-energy and vertex corrections as well as the renor- 
malization factors contain double poles such as 1/e 2 , 
l/(e') 2 and 1/ee'. However, we find that terms of the 
form - In all cancel out after we carry out expansions 
such as 



(-^) e = l + eln^ + 0(e 2 ). 

— LO —LO 



(A28) 



This is in accordance with the requirement of renormaliz- 
ability (2^^,^l|,^|. In addition, we find that the double 
poles in the renormalization factors cancel out in the beta 
functions. This is also a consistency check for the mini- 
mal subtraction procedure, and is in accordance with the 
requirement that beta functions are analytic. 



APPENDIX B: CALCULATIONS OF THE LOCAL 
SPIN SUSCEPTIBILITY 

1. The local susceptibility to order e 2 

The local spin correlation function is defined in 
Eq. (0). We find it convenient to carry out our calcu- 
lation in r space. (The self-energy and vertex correction 
diagrams in Appendix A were calculated in lo space.) So 
we take the expressions for the boson and conduction 
electron propagators as 



g°w = n„ (It 



sin 



(Bl) 







where Kq and Nq can be determined from Kq and A<o by 
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FIG. 12. Diagrams for the local spin susceptibility to the 



g order. 



Kl = K 2 T{2 - e) 
2 r(2-2e') 



m = n, 



(B2) 



°r 2 (l-e')' 

where T(x) is the Euler Gamma function. 

We now calculate the bare spin susceptibility in terms 
of the bare Hamiltonian (instead of the renormalized 
Hamiltonian used in Appendix A). The zero-th order 
result is given by 

X { b\t)=-\ fdr G%T)G%-T) = \e 



Xi 0) (r)= lim - X b(t) 

A— >-qo Z 



1 

4' 



(B3) 



where G°(r) is the bare /—electron propagator. This 
result is simply 5(5 + l)/3 with 5=1/2. 

To the order g 2 and J 2 , the diagrams have been shown 
in Fig. 8 of Ref. pTsj ] . Using the dimensional regulariza- 
tion adopted here, we find 



1 2(A g B ) V 
4 e 



l + Ce + \(C 2 +*-)e 2 



1 (N J B )W 7T 2 

'4 2e' U 6 6 h 



(B4) 



where C = 0.577215 • •• is the Euler Gamma constant. 

To g 4 order, the diagrams are shown in Fig. [l2[ (Struc- 
turally, they are similar to those considered in ref. [p4j| .) 
The result is 

^2 * 



1/4 1 + 8C 



2 + 2C + 8C^ 



e 



7T' 
~3 



x(A o5B )V 



(B5) 



We now express xb in terms of the renormalized cou- 
pling constants, g and J, by using Eqs. d9|jLl"|) in combi- 
nation with our results for the renormalization factors, 



Eq. (U). We can then determine the renormalization 
factor Z, defined in Eq. (|l8|), by demanding that all the 
pole terms in \b are removed by those of Z~ x so that 
the renormalized susceptibility, \ = % Xb> is regular. 
The result is 



x _ 2(K g) 2 (N Q J) 2 
Z - l + ^~ + ^2l^- 

| A(K g) 4 (K g) 4 



(B6) 



Eq. @ then follows. 

The renormalized susceptibility, \ = Z~ 1 xb, is equal 

to 

X (t) = X (0) [AO;, J, e) - J, e) ln^r) + 0(ln 2 (,ir))] , 

(B7) 

where, to the e 2 order, 

A(g, J, e) = 1 — 2(A 0ff ) 2 (C + i(C 2 + ^)e) 



+ (A o5 ) 4 (2 + 2C + 4C 2 - — ) 



(B8) 



J, e) = 2(A o3 ) 2 (l + Ce) + (N J) 2 - 2(A o5 ) 4 (l + AC). 

At the critical point("C" in Fig. [j]), the spin correlation 
has the form 



x(T)=X {a) A(g*,J*,e)^, 



(B9) 



where the exponent r\ can either be calculated from 
the renormalization factor Z in Eq. (|2l]), or by 
B(g*,J*,e)/A(g*,J*,e), which also yields r) = e. The 
renormalized amplitude factor A(g*, J*, e) is given by 

A{g\ J*, e) = 1 - Ce + Q + ic + ^C 2 - ^tt 2 ) e 2 . 

(BIO) 

For large values of e, it would require an appropriate 
resummation such as Pade approximation. 

For the local moment fixed point, the local spin corre- 
lation has the form 

X (r) = X (0) [A L (e) - B L (e) ln( M r) + 0(ln 2 (/ir))] , 
A L (e) = A(g,J,e)\ g=g * Li j = j* ==0 

B L {e) = B(g,J,e)\ g =g* L j=j l= o 
= e~Ce 2 . 

Again, B^e) / A^e) = e to the e 2 order. 



(Bll) 
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FIG. 13. The full local spin correlation function. Each dou- 
ble line represents the full /—electron propagator, while a 
shaded area denotes the full vertex. 



z n 



1 f(K g ± ) 2 {K g z f 



< \ 2 4 
1 fH Kpg^ (K g ± ) 2 (K o9z ) 2 (K o9z f 
8 8 32 



+ 



1 f(K o9± ) 4 (K og± ) 2 (K g z f 



e V 8 4 
1 f(N J^) 2 (N a J z ) 2 



e' V 8 



1G 



(C3) 



2. The local susceptibility to all orders in e 

We now discuss the contributions to x> to all orders of 
the perturbation theory. We can group all the diagrams 
in a manner illustrated in Fig. [l3|. Here, each double line 
represents the full /—electron propagator, while a shaded 
area denotes the full vertex. The renormalization factor 
for the spin, Z , is the product of the renormalization fac- 
tors for the full f-electron propagators and those for the 
vertices. The renormalization factor for a full f-electron 
propagator is equal to Zf. By inspecting the diagrams 
for each vertex, it is straightforward to see that to each 
order, this vertex is identical to that of the g— vertex. As 
a result, each vertex contributes a factor Z~ x . Taking 
these two things together, we end up with Eq. (p4[ ) , to 
all orders. As discussed in Section IV B, Eq. (|24D then 
leads to Eqs. @§7|). 



APPENDIX C: RG ANALYSIS FOR 
ANISOTROPIC CASES 

In the anisotropic case, the diagrams are topologically 
the same as those in the isotropic case, except for differ- 
ent coupling vertices. In the following, we will simply list 
the results. 

The renormalization factors are given as 

±{(K g±_) 4 , 3(K o9 ±) 2 (K o9z ) 2 (K o9z y 

2 



e° \ 8 8 
1 ( {IUg^ (K q9± ) 2 (K o9z ) 
e \ 8 4 

1 f(N J±) 2 , (N J Z ) 2 



32 



e' V 8 



1G 



(CI) 



V = 1 



(K g z ) 
4e 

J_ f (K Q gi 
e 2 \ 

(Ko9±) 4 , 



) 2 (K o9z ) 2 , {K o9z ) A 



4 

(N J Z ) 2 
16e' ' 



32 



(C2) 



Z J± = 1 - 



N J Z 



(K Q9z ) 2 , (N J Z ) 2 ( (N J ± ) 2 + (N J Z ) 2 



Zr = 1 



4e 16e' 2e' 2 

(K o9± ) 2 (N J z ) (K o9z ) 2 (N J z ) 

e'(e + e') 4ee' 
1_ f (K a9± ) 2 (K 09z ) 2 (K o9z ) A 
e 2 \ 4 32 

(K o91 _y 

8e 

N Q (Jl/J z ) , 1 f(K o9l _) 2 (K a9z f 



(C4) 



+ ^_f(N J ± ) 2 (N J Z ) 2 \ , (A J ± ) 2 



e' \ 8 16 
(K o9z ) 2 (N J z ) 

e'(e + e') 
(2(A g ± ) 2 - (K Q9z ) 2 ) (N Q jj/J z ) 
4ee' 

1_ / 3(A 0g± ) 4 _ (Kq 9 ^ 2 (K q9z ) 2 m (K 09z y- 
e 2 V 8 

1 f( A 'o.g±) 4 (A'n//_)-(Ao/y-.) : 
c 



32 



(C5) 



Eqs. (|C1|)- ( |C5| ) reduce to the isotropic results, Eq. (|12j), 
when we set <7 2 = <7j_, J 2 = J±. Similar to the isotropic 
case, we obtain the beta functions for each coupling con- 
stant from these renormalization factors, which are given 
in Eqs. @- @. We can then set g z = (gu = 0) to 
discuss the xy (Ising) case. 

We now turn to the calculation of the local spin suscep- 
tibility, focusing on the xy case only. Like in the isotropic 
case, we use the bare Hamiltonian. The contributions 
from the perturbative diagrams to the g 2 , J 2 , and g A 
orders are 

+_ (2) l (K o9 f) 2 + (K o9 ?) 2 e 

X B V) - 4 " T 

x ^ + Ce+ 1 -(C 2 + ^)e 2 

_{N Jlf + [ N jff , ^ 

We' 6 ; v ' 
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+-(4) M _ l (K gfY + (K,gf) 2 (K^) 2 2e 
Xb \ t > — 7 2 

e 3 



(There are no contributions at the J or g 2 J order.) Fol- 
lowing the same RG procedure as in Appendix B, we can 
calculate the renormalization factor Z: 



1 



1 



((K g ± ) 2 + (K g z ) 2 



i ({Kmy_ (WW 

e ^ 2 2 

1 f(N J±) 2 , (A J Z ) 2 



(C7) 



e' V 4 4 
and the anomalous dimension at the critical point: 
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l3i=9* ,Ji=J* 



dlogZ 
d\og[i 

(K gl) 2 + (K g* z ) 2 - (Kogl) 4 

(k *\ 2 ik *^,(n ji) 2 , (iv j;) 2 



(C8) 



2 2 

Setting g z = (for the xy case) , the exponent becomes 
(v *^2 ( V *\4 . (iVo J^) 2 + (7V J*) 2 

v=[Kog±) -{Kog±) + ^ ■ (C9) 

The renormalized amplitude A(g, J, e) is given by 

A(g, J, e) = 1 - (^o.9±) 2 (C + ^(C 2 + ^)e) 

+ (^ .9±) 4 (1 + C + 2C 2 -^) (CIO) 

6 

At the critical point, taking account of Eq. (f37j), we 
find 

»?= e, 

A c (e)= 1 - Ce + (l + + ^C 2 - i. 2 ) £ 2 . (Cll) 

Near the bosonic stable fixed point where = J* = 0, 
which is given by Eq. j39|), the results are 

riL = (K gl) 2 - (Kogl) 4 
= £) 

^ L ( e ) = l-C e +(l + ^ 2 - V) £ 2 . (C12) 



given in ref. The resulting action, describing a one- 
dimensional long-ranged statistical-mechanical model, is 
similar to Eq. (7) of ref. H] and, in our notation, takes 
the following form: 



S(r 2n , ri) - -2n ln( W ) + ft (n+i - n)/ 



TO 



J2(-l) l+m P*i Hr m - n)/r + K(r m - n)] (Dl) 



where [T% n , ...,ri], for n = 1,2,..., labels a sequence of 
spin flips (kinks) along the imaginary time axis, and yj 
and Kj are defined in Eq. ©. The last term, K(r), 
originates from the g z — coupling and takes the form: 



K(r) = Kg [(r/r o r-l}/e 



(D2) 



where K g is also specified by Eq. (|40|). The RG equations 
for this kink-gas problem has already been derived in 
ref. EM and are given as follows: 



[3( Kj ) 
p( Kg ) 
13(h) 



- Vj (l - Kj - Kg/2) 

\k jV 2 

-K g (e-4y 2 ) 
-Hi - 2y 2 ) 



(D3) 



This procedure is valid for arbitrary values of the stiffness 
constants Kj and K g , provided the fugacity yj is small. 

An unstable fixed point still exists. To the linear order 
in e, it is located at 



(y 



*\2_ 



6/4 

2 



K*=0 



(D4) 



From Eq. (D3), we can easily see that K g is in fact 
irrelevant around this fixed point. We can then determine 
the critical exponents by staying within the yj — K g plane: 
the projection of the RG flow on the yj — K g plane is given 
in Fig. ||. Within the e— expansion, then, the critical 
exponents of the unstable fixed point become identical to 
their counterparts for a similar fixed point of the classical 
ferromagnetic Ising chain with an long-range interaction 
that decays in distance in terms of a power-law exponent 
2 — e JO]]. The critical exponent 77, for instance, can be 
strai ghtf orwardly calculated by combing Eq. ( |D3j ) and 
Eq. ( D4 ) . The result is r\ = e, the same value as it takes at 
the critical point of the corresponding long-ranged Ising 
chain fltl. 
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